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ABSTRACT 

We present the first results of a sub-millimeter survey of distant clusters using the new Sub-mm 
Common-User Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope. We have mapped 
fields in two massive, concentrated clusters, A370 at z — 0.37 and CI 2244—02 at z = 0.33, at wavelengths 
of 450 and 850 /im. The resulting continuum maps cover a total area of about lOarcmin 2 to la noise 
levels less than 14 and 2mJybeam _1 at the two wavelengths, 2-3 orders of magnitude deeper than 
was previously possible. We have concentrated on lensing clusters to exploit the amplification of all 
background sources by the cluster, improving the sensitivity by a factor of 1.3-2 as compared with a 
blank-field survey. A cumulative source surface density of (2.4 ± 1.0) x 10 3 degree -2 is found to a 50% 
completeness limit of ~ 4 mJy at 850 /im. The sub-mm spectral properties of these sources indicate 
that the majority lie at high redshift, z > 1. Without correcting for lens amplification, our observations 
limit the blank-field counts at this depth. The surface density is 3 orders of magnitude greater than 
the expectation of a non-evolving model using the local IRAS 60 (im luminosity function. The observed 
source counts thus require a substantial increase in the number density of strongly star-forming galaxies 
in the high-redshift Universe and suggest that optical surveys may have substantial underestimated the 
star formation density in the distant Universe. Deeper sub-mm surveys with SCUBA should detect large 
numbers of star- forming galaxies at high redshift, and so provide strong constraints on the formation of 
normal galaxies. 

Subject headings: cosmology: observations — cosmology: early universe — galaxies: evolution — galax- 
ies: formation — gravitational lensing — radio continuum: galaxies 



1. INTRODUCTION 

Surveys of the local Universe have shown that a third 
of the total bolometric luminosity is emitted at sub-mm 
and far-infrared wavelengths as a result of reprocessing of 
star-light by dust (Soifer & Neugebauer 1991). Moreover, 
some of the most vigorous star-forming galaxies in the local 
Universe are also those in which the effects of dust obscu- 
ration are most significant. While there have been striking 
advances in the identification of 'normal' galaxies at high 
redshift {z ~ 2-4) using Lyman-dropout techniques (Stei- 
del et al. 1996), such approaches are insensitive to highly 
obscured star-forming galaxies at these epochs. The pres- 
ence of at least modest amounts of dust in distant proto- 
galaxies, especially forming spheroids, is expected given 
the highly metal-enriched ISM which must be present dur- 
ing their formation (e.g. Mazzei & de Zotti 1996). Thus 
techniques sensititive to obscured distant galaxies should 
be used to investigate the formation of these populations 
in the early Universe. 

Sensitive sub-mm observations present the first opportu- 
nity to detect dusty star-forming galaxies at high redshift. 
At wavelengths around 100 /im, the bulk of the luminosity 
of normal, star-forming galaxies is reprocessed star-light 
from dust and so observations in the sub-mm band can 
provide robust estimates of both the dust mass and total 
star formation rate in a galaxy. Furthermore, the negative 
if-correction at wavelengths A ^ 400 fim means that sub- 
mm observations select star-forming galaxies at z <; 1 in 
an almost distance-independent manner, providing an effi- 
cient method for the detection of star-forming galaxies at 
very large redshifts, z < 10, and hence studying galaxy 



evolution (Blain & Longair 1993, 1996 — BL96; Blain 
1997; Eales & Edmunds 1997; Franceschini et al. 1997; 
Guiderdoni et al. 1997). 

Most published sub-mm studies of distant galaxies have 
targeted atypical galaxies (e.g. radio-loud galaxies, Ivison 
et al. 1998). We report here the first deep sub-mm survey 
to probe the nature of normal galaxies at moderate and 
high redshift, z > 0.5. In this study we have attempted to 
maximise the available sample of distant galaxies by con- 
centrating on fields in moderate-redshift clusters. While 
the dominant spheroidal populations of these clusters are 
expected to be quiescent in the sub-mm band, the in-fall 
of field galaxies associated with the growth of the clusters 
(Smail et al. 1997) means that these fields may contain 
over-densities of moderate-redshift star-forming galaxies, 
as compared with 'blank' field surveys. 

The main attraction of the clusters observed here, how- 
ever, is that they are strong gravitational lenses, magni- 
fying any source lying behind them (Blain 1997). Given 
the expected steep rise in the sub-mm counts (BL96), am- 
plification bias could increase the source counts in these 
fields by a substantial factor, with a maximum predicted 
surface density of about 10 sources per SCUBA field down 
to 1 mJy at 850 /im (Blain 1997). Moreover, by targeting 
those clusters that contain giant arcs, images of distant 
field galaxies magnified by factors of 10-20, we can also 
obtain otherwise unachievable sensitivity (< 0.1 mJy at 
850 /im) on the dust properties of a few serendipitously- 
positioned normal galaxies at high redshift. The angular 
scales of the region where these highly magnified high- 
redshift galaxies are found is also well-matched to the 
SCUBA field-of-view. 
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Tabic 1 

SCUBA observations of A370 and CI 2244-02 



Target 


R.A. 


Dec. 


A 


Exposure time 


Area 


Flux density 




(.T2000) 


(J2000) 


(/"») 


(ks) 


(arcmin 2 ) 


la (mjy) 


CI 2244-02 


22 47 11.9 


-02 05 38 


450 


23.0 


4.00 


14.0 








850 


23.0 


5.40 


1.9 


A370 


02 39 53.0 


-01 35 06 


450 


25.7 


4.35 


13.3 








850 


25.7 


5.40 


1.8 



In the following sections we give details of the observa- 
tions and their reduction, and discuss the results within 
the framework of current theoretical models of galaxy for- 
mation and evolution. We adopt H = 50 km s _1 Mpc -1 
and q = 0.5. 

2. OBSERVATIONS AND REDUCTION 

These data were obtained using SCUBA (Cunningham 
& Gear 1994) on the James Clerk Maxwell Telescope 
(JCMT) 5 . SCUBA contains a number of detectors and de- 
tector arrays cooled to 0.1 K and covering the atmospheric 
windows from 350 /xm to 2000 /im. In our survey, we oper- 
ated the 91 element Short-wave (SW) array at 450 /im and 
the 37 element Long- wave (LW) array at 850 /xm, giving 
half-power beam widths of 7.5 and 14.7arcsec respectively. 
Both arrays have a 2.3 arcmin instantaneous field-of-view 
and the design of the optics ensures that, with a suit- 
able jiggle pattern for the secondary mirror,, fully sampled 
maps can be obtained simultaneously at 450 and 850 xtm. 
The multiplexing and high efficiency of the arrays means 
that SCUBA offers a gain in mapping speed of a factor of 
about 300 as compared with previous detectors. 

The observations employed a 64-point jiggle pattern, 
fully sampling both arrays over a period of 128 s. The 
pattern was subdivided so that the target position could 
be switched between the signal and reference every 32 s: a 
repeating signal-reference-reference-signal scheme, with 
sixteen 1 s jiggles in each. Whilst jiggling, the secondary 
was chopped at 6.944Hz by 60arcsec in azimuth. The 
pointing stability was checked every hour and regular sky- 
dips were performed to measure the atmospheric opacity. 
The rms pointing errors were below 2 arcsec, while atmo- 
spheric zenith opacities at 450 and 850 /xm were very sta- 
ble during the course of each night, the night to night 
variations being in the range 0.98-2.2 and 0.18-0.37 re- 
spectively. 

The dedicated SCUBA data reduction software (SURF, 
Jenness 1997) was used to reduce the observations. The 
reduction consisted of subtracting the reference from the 
signal after carefully rejecting spikes and data from noisy 
bolometers. Six quiet bolometers at the edge of each array 
were used to compensate for spatially-correlated sky emis- 
sion. This reduced the effective noise-equivalent flux den- 
sity from 100-350 to gOmJyHz- 1 / 2 at 850 /xm (c.f. Ivison 
et al. 1998). The resulting maps were fiatfielded, corrected 
for atmospheric attenuation, and calibrated using nightly 
beam maps of Uranus. The calibrated maps from each 
night were coadded and then linearly interpolated onto 
an astromctric grid using an approximately Nyquist sam- 
pling, of 2 and 4 arcsec pixel -1 at 450 and 850 /xm respec- 
tively, to produce the maps presented in Fig. 1 (Plate 1). 

The final on-source integration times are listed in Ta- 
ble 1, along with field positions and sensitivity limits. The 
beams have moderate error lobes, especially at 450 /xm, 



and so we correct our aperture measurements for flux 
outside the aperture using the measured values from our 
calibration sources. Even without including a factor to 
account for the lensing amplification, the data shown in 
Fig. 1 are the deepest sub- mm maps ever published, and 
illustrate the cosmetically clean and flat maps achievable 
with SCUBA in long integrations. 



3. ANALYSIS AND RESULTS 

Source catalogs from our fields were constructed using 
the Sextractor package (Bcrtin & Arnouts 1996). The de- 
tection algorithm requires that the surface brightness in 4 
contiguous pixels exceeds a threshold (chosen as ~ la of 
the sky noise, Table 2), after subtracting a smooth back- 
ground signal and convolving the map with a 4 x 4 pixel 
top-hat filter. The numbers of objects (N) detected in each 
field are given in Table 2, indicating that our observations 
are far from being limited by confusion: there are ~ 60 
beams per source. 

To assess the contribution of noise to our catalogs we re- 
ran the detection algorithm on the negative fluctuations 
in the map. This gave a simple estimate of the number 
of false-positive detections that may arise from the noise, 
assuming that the noise properties of the map are Gaus- 
sian. We estimate that there are no false detections in 
our catalogs (N_„ e , Table 2), and so all the detections are 
real. The presence of the brightest source in the refer- 
ence beams (60 arcsec to the East and West in the 850 /xm 
map of A370) was disregarded. This detection does con- 
firm, however, the reality of positive features at this faint 
level, while the absence of any other negative detections 
limits the number of luminous sources which can lie in the 
regions covered by the reference beams. 

Secondly, to determine the completeness of our sam- 
ple we added a faint source to the maps repeatedly, re-ran 
our detection algorithm and estimated the efficiency of de- 
tecting this source as a function of its flux density. This 
provides a reliable estimate of the visibility of a faint com- 
pact source in the maps. The template source was a scaled 
version of our calibration source, Uranus. The estimated 
80% and 50% completeness limits of the catalogs derived 
from these simulations are listed in Table 2 as S 80 % and 
$50%- The incompleteness limits are relatively bright for 
the 450 /xm maps because a large proportion of the flux 
density (about 40%) is found in the low surface brightness 
wings of the beam. The simulations also indicate that 
the measured 850 /xm flux densities are unbiased and are 
typically accurate to 10% at 25mJy and 30% at 4mJy. 

We discuss the detailed properties of the sources in an- 
other paper, where we also place limits on the dust masses 
of the numerous strongly-lensed distant galaxies covered 
by our maps. However, we note that, based on their weak 
or non-detection at 450 /xm, all of the 850 /xm sources in 

5 The JCMT is operated by the Observatories on behalf of the UK Particle Physics and Astronomy Research Council, the Netherlands 
Organization for Scientific Research and the Canadian National Research Council. 
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Table 2 
sub-mm source counts 



Target A Threshold Min. Area N N_„ c S 80 % S 50 % 

(/^m) (mjy beam -1 ) (arcscc 2 ) (mjy) (mjy) 

CI 2244-02 450 16.0 16 135 63 

850 1.9 64 2 5.5 3.9 

A370 450 16.0 16 1 140 84 

850 1.8 64 4 5.4 3.9 



our sample appear to have the sub- mm spectral character- 
istics of distant (z <; 1) star-forming galaxies and are thus 
unlikely to be associated with the clusters. 




1 2 5 10 30 



Limiting Flux Density / mjy 
Figure 2 Models of the integral number counts of sources at 850 /^m in 
a parametric model of galaxy evolution (BL96) and from a simple model 
based on the limits on strongly star-forming systems in optical surveys 
of distant galaxies. The observations arc represented by filled circles 
with Poisson errors, note that the errors arc not independent on the var- 
ious points. The observations have been corrected for the effects of lens 
amplification using simple lens models, but not corrected for incomplete- 
ness. The solid curves represent, in order of increasing predicted counts, 
models that include: no evolution; (1 + z) 3 evolution with -s ma x — 2 
and zq — 5 (Model A); and (1 + z) 3 evolution with z max — 2.6 and 
zq — 7 (Model B). The dashed lines represent models where we fill the 
Universe across z —0-10 with a constant density (0.6 X 10 4 Mpc~ 3 ) of 
star-forming galaxies with fixed star formation rates. In order of increas- 
ing predicted counts the dashed lines represent star formation rates for 

the population of: M— 20, 50 and 150 MQyr -1 , where we have assumed 
a dust temperature of T — 60 K. Clearly only models including high den- 
sities of strongly star-forming galaxies are compatible with the observed 
surface density of sources. 

Converting the observed number of sources into a sur- 
face density and correcting for incompleteness, we deter- 
mine a cumulative number density across our two fields of 
(2.4 ± 1.0) x 10 3 degree -2 down to a 50% completeness 
limit of 4 mjy at 850 /xm (all errors include only Poisson 
contributions). At 450 /im, the single source we detect 
places only a very weak limit on the likely surface density 
of the order of 1 x 10 3 degree -2 brighter than 80mJy. 

We now estimate the likely lens amplification factors, 
and so place tighter limits on the typical blank-field counts. 
Because the distances to the detected sources are un- 
known, this estimate will, by necessity, be crude and so 
we have not attempted a detailed analysis. The cluster 
potentials are modelled as isothermal spheres with masses 
and centers determined from the redshifts and observed 
shapes of the giant arc in each cluster (Kneib et al. 1993; 
Smail et al. 1996). In these models the mean amplification 
factors for typical background sources (z > 1) are about 
2 and 1.3 in the regions covered by our maps of A370 
and CI 2244— 02 respectively, while the observed area of 
the maps (5.4 arcmin 2 at 850 /im) corresponds to 1.8 and 
4.0arcmin 2 in the respective source planes. Correcting 
the flux densities of our sources to take account of the 
probable lens amplifications, but not correcting for incom- 
pleteness, we predict the source counts presented in Fig. 2. 



These indicate integrated number densities in blank fields 
of (2.5 ± 1.4) x 10 3 (similar to the observed value due 
to the source counts having a form close to S 0A ) and 
(3.6 ± 1.6) x 10 3 degree -2 to flux limits of 4 and 3mJy 
respectively at 850 /im. 

From the flux densities associated with the re- 
solved sources in the fields we calculate lower limits to 
the background radiation intensities of 2.6 and 2.4 x 
10~ 10 Wm -2 sr -1 at wavelengths of 450 /mi and 850 /an 
respectively, averaged over both fields. By extrapolating 
the 850 /mi counts using our best-fit model (Fig. 2 and 
§4) to faint flux densities, we estimate total intensities of 
extragalactic background radiation from discrete sources 
of about 26-28 and 4.4-6.7 x 10 -10 W m -2 sr" 1 at wave- 
lengths of 450 and 850 /im respectively. These background 
radiation intensities are broadly consistent with the tenta- 
tive detection of an isotropic component of the background 
radiation in the sub-mm by Puget et al. (1996), who in- 
ferred vl v = (23±20) and (2.7±2.0) x 10 -10 W m -2 sr -1 at 
wavelengths of 450 and 850 /mi respectively. If we assume 
that all of the background radiation intensity we infer is 
due to the formation of massive stars, then we expect that 
a density parameter of heavy elements of £ 6 x 10 -4 will 
have accumulated in the Universe by the present epoch. 
This density corresponds to about 1.1% of the density pa- 
rameter in baryons if fib = 0.05, and so it is fully consis- 
tent with present limits. We reiterate, however, that these 
are tentative estimates, the accuracy of which depends on 
the models assumed for both the lens and the form of the 
counts of distant galaxies. 

4. DISCUSSION 

Due to the negative if-corrections expected for distant 
galaxies, sub-mm observations provide a good estimate of 
the volume density of luminous star-forming galaxies at 
z > 1. In the absence of redshifts for all the sources, the 
evolution of this population can be understood by com- 
paring parameterised models to the counts (BL96). The 
BL96 models are based on the 60 /im luminosity function 
of IRAS galaxies (Saunders et al. 1990) and assume that 
the luminosities evolve as (1 + z) 3 out to a rcdshift, z max , 
and then maintain the enhanced luminosity out to a cutoff 
rcdshift, zq. The form of this evolution is motivated by the 
observations of similar behaviour in both the radio galaxy 
and QSO number counts (Dunlop & Peacock 1991) as well 
as the luminosity density of field galaxies at z < 1 (Lilly 
et al. 1996). BL96 also give predicted counts for a non- 
evolving model using the same luminosity function. The 
adopted parameters for the models described in that paper 
give source counts which roughly span the range predicted 
by other similar works (e.g. Guidcrdoni et al. 1997). In 
Fig. 2, we plot both the no-evolution case and two other 
parametric models based on BL96: model A - Model 2 in 
BL96 - with values of z max = 2 and zq = 5; and model B, 
has z max ~ 2.6 and z = 7, although most combinations of 
£max — 2.2-2.9 and z k, 5 give comparable results. Model 
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B was used to estimate the total background radiation in- 
tensity discussed in §3. 

From Fig. 2 it can be seen that the no evolution predic- 
tions fall short by 2-3 orders of magnitude of the obser- 
vations. Thus, this first analysis of a deep sub-mm survey 
indicates that the number density of strongly star-forming 
galaxies, and so the mean star formation rate in the dis- 
tant Universe is considerably larger than that seen locally. 
To estimate the extent of this evolution we assume that 
all the detected sources lie beyond the clusters. We then 
require strong evolution, of the form given in model B, 
out to z > 2 to fit the 850 /xm counts. For consistency, we 
check the predictions from model B for the observed counts 
at 450 fim; 0.7 sources are expected in the two fields, in 
agreement with the single detection. 

We conclude from the 850 ^m counts that the integrated 
star formation rate in the Universe, as traced by the num- 
ber density of the most luminous sources, must continue 
to rise out to beyond z > 1 (z > 2 in our parametric 
models), extending the trend observed at z < 1 (Lilly 
et al. 1996). For sources at z ~& 1, the typical luminos- 
ity we infer is Lfir ~ 0.5-1.0 x 10 13 L Q , with a star for- 
mation rate of 100-300M©yr _1 . Using the observed 
surface density of these objects and assuming a constant 
space density of sources between z = 1-5, we would pre- 
dict a number density of strongly star-forming galaxies of: 
N{MZ ISOMoyr- 1 ) - 1.2 x 10~ 4 Mpc~ 3 , at z £ 1. 

Limits on the number density of strongly star-forming 
galaxies at z ~ 2-3.5 have recently been published by 
Madau et al. (1996 — M96) on the basis of Lyman- 
dropout surveys. Their limit is N(M> 20M Q yr _1 ) < 
0.6 x 10~ 4 Mpc~ 3 . We plot in Fig. 2 three models us- 
ing this number density of sources to uniformly populate 
the volume from z = 0-10, but allowing the corresponding 
star formation rate to vary: M= 20, 50 and 150MQyr _1 . 
These values represent that originally given by M96, the 
M96 limit corrected for dust extinction as suggested by 
Pettini et al. (1997), and a mean star formation rate closer 
to that needed to fit our observations. A galaxy popu- 
lation consistent with the limits from Madau et al. pre- 
dicts a source density 3 orders of magnitude lower than 
that observed (Fig. 2). Even if the modest dust extinc- 
tion proposed by Pettini et al. is included, we still under- 
predict the observed surface densities (unless the dust in 
this population is very cold, T = 40 K, and they have very 
large dust masses). To match the observed surface den- 



sity of 850 /im sources we must significantly increase the 
mean star formation, cither by further increasing the star 
formation rate associated with the optically-selected sam- 
ples (c.f. Meurer et al. 1997) or by introducing a popula- 
tion of strongly star-forming, but highly obscured, distant 
galaxies missed in these samples. Deeper sub-mm surveys 
(BL96; Pearson & Rowan-Robinson 1996) are necessary to 
differentiate between these possibilities and hence provide 
an unbiased view of star formation in the distant Universe. 

5. CONCLUSIONS 

• We have presented the first sub-mm survey of the distant 
Universe, deep enough that we should detect the evolving 
galaxies predicted by current theoretical models, while at 
the same time covering a sufficiently large area to be sta- 
tistically reliable. We derive cumulative source counts of 
(2.4 ± 1.0) x 10 3 degree -2 down to 4mJy at 850 ^m. 

• The surface density of faint sources in the sub-mm far 
exceeds a simple non-evolving model using the locally ob- 
served 60 fj,rn galaxy luminosity function. Thus our obser- 
vations require a substantial increase in the number den- 
sity of strongly star- forming galaxies at z > 1. 

• Comparison of our observations with the predictions of 
simple parametric models implies that the luminosity den- 
sity of the brightest sub-mm sources continues to increase 
out to z > 1. Models based upon the claimed properties 
of star-forming galaxies from optically-selected samples of 
distant galaxies (Madau et al. 1996) significantly underes- 
timate the observed surface density of sub-mm sources. 
We suggest that such samples may be missing a consid- 
erable proportion of the star formation in dust-obscured 
galaxies at high redshift. We conclude that question of 
the evolution of the star formation density in the distant 
Universe and hence the epoch of galaxy formation is still 
very much open. 
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Figure 1 (Plate 1). The 450 and 850 fim maps of the two fields: a) A370, 850 /im; b) CI 2244—02, 850 /im; c) A370, 450 ^m; d) CI 2244-02, 450 fim. 
The maps are smoothed to the instrumental resolution at each wavelength and have had their boundaries apodized to remove regions on the outskirts 
with low effective exposure times resulting from the rotation of the field-of-view during our long exposures. They are displayed as a grayscale from —4a 
to 4a, the contours are positive and show 3, 4, 5, 10, 15 a for each field. The major tick marks are 20 arcscc in all panels. 



